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Abstract 

We consider the quantum mechanical hamiltonian of two, space 
indexed, hermitean matrices. By introducing matrix valued polar co- 
ordinates, we obtain the form of the laplacian acting on invariant 
states. For potentials depending only on the eigenvalues of the radial 
matrix, we establish that the radially invariant sector is equivalent to 
a system of non interacting 2 + 1 dimensional fermions, and obtain its 
density description. For a larger number of matrices, the presence of 
a repulsive radial inter-eigenvalue potential is identified. 
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1 Introduction 



The study of multi-matrix models^, and particularly their large N limit pQ, 
is of great interest. It is well known, for instance, that the large N limit of 
their description of DO branes [2] has been conjectured to provide a definition 
of M theory [3]. In the context of the AdS/CFT duality [I], [5], [6], due to 
supersymmetry and conformal invariance, correlators of supergravity and 
1/2 BPS states reduce to calculation of free matrix model overlaps [7J, [8] or 
consideration of related matrix hamiltonians For stringy states, in the 
context of the BMN limit [10] and M = 4 SYM, similar considerations apply 
HU, H3, [33]. A plane- wave matrix theory [H] is related to the Af = 4 SYM 
dilatation operator [15] . Recently, multi-matrix, multi-trace operators with 
diagonal free two point functions have been identified [16], [17] . In earlier 
works [18], [19], [20] it has been argued that QCD can be reduced to a finite 
number of matrices with quenched momenta. 

As such, the study of the quantum mechanics of two hermitean matrices, 
or equivalently of a single complex matrix, is of particular interest, as this 
system already possesses many of the complexities of multi-matrix models. 

Previous related studies of two hermitean matrices have typically treated 
the two matrices asymmetrically. For instance, in the approach first de- 
veloped in [21J, one of the matrices, which generates the large N planar 
background, is treated in a coordinate representation, and the other in a 
creation/anihilation basis. This was done in the context of 1/2 BPS states 
and a dual free harmonic Hamiltonian [8] , [9] with the matrix generating the 
background being the holomorphic component of a complex matrix. A pre- 
cise phase space identification between the collective density description of 
the dynamics of this matrix [22] , and the droplet description of the LLM [23] 
metric is obtained |21j. The generalization of this approach to include gyM 
interactions was developed in [24J. By considering the planar background 
generated by of one of the two Hermitean matrices, further properties of 
the spectrum were established in [25]. Non supersymmetric, gyM dependent 
backgrounds have been considered in [26], but still within an aproach where 
the two matrices are treated asymmetrically. 

In Sections 2 — 5 of this communication, we consider a more symmetrical 
approach to the quantum mechanics of two hermitean matrices X\ and X2, 

1 By matrix models we mean integrals over matrices or the quantum mechanics of matrix 
valued degrees of freedom 
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with 1,2 spatial indices (or complex matrix Z — X 1 + iX 2 ), by introducing 
matrix valued polar coordinates. In particular, we study the kinetic energy 
operator, or laplacian. In Section 2, after diagonalizing the "radial" matrix, 
two sets of variables are identified in terms of which it is possible to invert the 
metric and obtain its determinant. For both sets, there is no mixing between 
the radial eigenvalue differentials and the angular differentials. Both metrics 
depend on the eigenvalues of the radial matrix. The corresponding laplacians 
are derived. In Section 3, we obtain the unique laplacian that acts on gauge 
invariant states, expressed in terms of the eigenvalues of the radial matrix and 
of a single "angular" unitary matrix. A constraint on the set of these states 
is identified. In Section 4, for potentials which depend only on the radial 
eigenvalues, we show that the radially symmetric sector is equivalent to a 
system of non- interacting 2 + 1 dimensional " s-state" fermions. In Section 5 
we consider the density description of this system, and obtain the additional 
potential resulting from the radial inter-eigenvalue repulsion. 

In Section 6, we consider an arbitrary number of complex matrices and 
identify in this case too, the presence of repulsive radial inter-eigenvalue 
interactions. We conclude with some comments in Section 7. 

2 Matrix polar coordinates 

We wish to consider the quantum mechanics of two NxN hermitean matrices 
Xi and X 2 with hamiltonian 

£ 1 ( d d d d 

~2 {diXjijdiXjjt + d{X 2 ) ij d{X 2 ) ji 

= ~V 2 + V(X 1 ,X 2 ) 

We introduce matrix valued polar coordinates 

X 1 + iX 2 = Z = RU , tf = U*R (2) 

with R hermitean and U unitary. Since R is hermitean, it can be diagonalized 
as R = V'rV, with r a diagonal matrix and V unitary, and we obtain 

Z = V*rVU = V^rW, W = VU 
Z^ = U^rV = W ] rV (3) 

3 



j+V(X 1 ,X 2 ) 

(1) 



We will refer to the parametrization in terms of (r, V, U) as parametrization 
(I), and the parametrization in terms of (r, V, W) as parametrization (II). The 
number of degrees of freedom is preserved, since these matrix coordinates are 
defined up to V — > DV, W — > DW with D a diagonal unitary matrix. 
Introducing the anti-hermitean, Lie-algebra valued differential matrices 

dX = VdUU^V ] , dS = dVV^ , dT = dWW* (4) 

we obtain: 



dZ = V ] {dr+[r,dS]+rdX)VU = V ] {dr + rdT-dSr)W 
dZ^ = UW(dr+[r,dS]-dXr)V = W\dr + rdS - dTr)V (5) 

The metric is defined from the infinitesimal length squared 

TrdtfdZ = Tr [{drf + [r, dS] [r, dS] - r 2 {dX) 2 + [r, dS] [r, dX}] 

= Tr [(dr) 2 -r 2 (dS) 2 -r 2 (dT) 2 + 2rdSrdT] (6) 

Starting with parametrization (I), we note that the infinitesimal length squared 
has a 'local'(or "pointwise") form in terms of the double index i.e.: 



TrdtfdZ = ^drl-Y^in-rjfdSijdSji 

i ij 

~ ^(ri-rtfidSydXii + dXvdSji} (7) 

ij 

~ lT,(r- + r^dX.jdXj, 

ij 

Recalling the antihermiticity of the differentials and writing ds 2 = gA,Bdx A dx B 
with dx A = (dri^X^dSij^^Xij^j^dS-^^^dX*^^), we obtain for 
9a,b- 



( 1 

r 2 













(ri-rj) 2 X( - --^ 2 



2\ r i- r 3> 



\iji-rj) 2 \(r 2 + r~ 










\ 




\ 








[Ti - Tj 



12 Wi-Tjf 



Wi-Tjf \{r 2 + r 2 ) J 
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It follows that 
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det g AB = l[r> (]J \{r 2 - r]A = {\\r^\ IR f 

i \i<j / i 



(8) 



where we have defined 



Aw^n^ 2 -^) 2 - ( 9 ) 

i<j 

The laplacian is obtained in the standard way, and it takes the form: 



v « = rT^A^7rfn^ A2 ^)^ ( 10 ) 

2(r 2 + r 2 ) d d 4 8 8 



x ^ 2 /_g d_ 8 8 



For parametrization (II), we note that the second expression for the infinites- 
imal length squared in equation ([6]) can be further diagonalized, for each i 
and j, by introducing: 

dY + = -^{dT + dS) dY~ = -^=(dT - dS) 
V 2 v 2 

Then 



TrdZ f rfZ = Tr[dr 2 + -[r, [r, - -{r, dy~}{r, dK"}] (11) 
Writing again ds 2 = TidZ^dZ = gA,Bdx A dx B ', we obtain 

det g A>B = X\2r 2 t {A) 

i 

and straightfowardly 
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v <-» - raM^M (12) 

2-^ 9^2 ay- fly*- + 



+ > 7 > 



3 Invariant states 

We are interested in the action of the above laplacians, or Hamiltonian, on 
invariant states, i.e., states obtained by tracing strings of Z's and Z^s : 

Tr(...Z np Z^ mp ...Z nq Z^ mq ...) 

These states depend only on the eigenvalues of the radial matrix R and 
on the unitary matrix 

Q = VUV^ = WV ] 
Starting with parametrization (II), we note that 



dQ = dTQ - QdS = -^={dY + Q - QdY + ) + -^{dY'Q + QdY~) 

v 2 y2 

from which, on this invariant subspace, 

9 <9Q ab d _ 1 (l) p (fl). 

where we have introduced the generators of left and right U(N) "rotations": 



" l:h <)(),,, " '-0Q 
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Substitution of these expressions into the laplacian (TT2|) yields 



v 2 = ^£^.A? ^ H)^-£^(4 , + 4 R, )(4 , + 4 R, ) (i4) 

MR i >%<-><% U'l i I' i 

_ 2(r ' 2 + ^ f £ (L) + ) - 4r ' rj f E [R) + E iR) E (L) )) 

Before arriving at the final expression, note from fill I) that dY^ is absent, 
and therefore we require: 

o = 5 | f = ^(4 1, -4' ,) ) ; 4 1) = 4 K) (15) 

The final form of the laplacian is then: 



^MR i r i OTi OTi ^ rf 

V ( 2 ( r i +r j\ E (L) E (L) E (R) _ 4r ' r i ( E (L) E {R) , E (R) E (Lh\ 

A -J- A \ V % J ' ^ i 3 ' * 



2 ( r i +r j) , F {L) F (L) p (R) F (R)- ] _ 4r » r i f jp(L) , p (Ji) 

^ ^ r 2 _ r 2)2 V «i "'" «7 J* / ( r 2 _ r 2 

In terms of parametrization (I), from 

dP = dQQ ] = dS + dX- QdSQ^ 
one obtains, when acting on invariant states, 



9 _ jp(L) 9 _ F (L) F (R) 



Substitution of these in ( flOl) results in the laplacian (1161) . with states hav- 
ing to satisfy the contstraint ( fl5l) . now resulting from the ciS^ mode. As 
expected, the form of the laplacian ffl6l) is unique. 

4 Radial fermionization 

In the case of the simple integral (d = system) over two hermitean matri- 
ces, it follows from the form of the '"volume"' element (JBJ) that an additional 
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term is added to the action expressing the the repulsive inter-egenvalue inter- 
action amongst the radial eigenvalues. By either integrating out the angular 
dependence, or by factoring it out if the action is only dependent on the ra- 
dial eigenvalues (as is the case for the free theory in an harmonic potential) 
[27], one straighforwardly obtains a BIPZ [28] type equation. 

For the hamiltonian system, it is a well known result [28J that the singlet 
sector of a N x N single hermitean matrix hamiltonian with a potential de- 
pending only on its eigenvalues is equivalent to a system of TV non-interacting 
fermions. This is a result of the anti-symmetry under the exchange of any 
two coordinates of the Van der Monde determinant 

A(x k ) = Y[( x i - x j) 

i<j 

We consider the case of a potential that depends only on the eigenvalues 
of the radial matrix U. An example would a potential of the form: 

V(X l ,X 2 ) = Trv(ZZ^) = Trv(Z^Z) 

with v(x) a polynomial. Then, on "s-states" (independent of the "angular" 
degrees of freedom Q), and letting pi = rf, 



"2 V = "2 Af^ ^ nfri [r ^ MR) frl = " A^) ^ Wi [Pl {P)) 

This kinetic energy operator acts on symmetric wavefunctions $. Defining 

^ = A$ , (18) 

its action on ^ takes the form: 
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The second term clearly vanishes, and the last two vanish due to the identity 



which generalizes the identity applicable to the single hermitean case. It is 
easily proven by choosing any three distinct eigenvalues. 

Therefore, the eigenvalue equation for the energies of the system takes 
the form, 



This is the "s-state" Schroedinger equation for N non-interacting 2 + 1 di- 
mensional non-relativistic fermions. 

5 Density description 

In this section we describe the key features of the density descripton of ra- 
dially symmetric fermions int terms of the density of radial eigenvalues. We 
use the collective field theory approach of Jevicki and Sakita [22] . 

The existence of such a descripton requires the identification of a suit- 
able set of invariant operators which close under "'joining'" and '"splitting"', 
equivalent to the closure of underlying Schwinger-Dyson equations. 

Remarkably, the following set can be identified as such: 



% ± k 




( d d , \ T / 1^1 9 d 

{- 2 ^W:"'W,' v{pi) ) * = [-2 E j^Wi 



(20) 




(21) 



One has 




(22) 
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from which the " 'joining" ' operator Q |22j takes the form 



dZl dZ 



-kk'TrZ^Ze l(k+k '^ z 



(23) 



j i 



dk f dk' 



e- ikx e- ik ' x Q kk , = d x d x , [x$(x)5(x - x')] 



2tt J 2tt 

For the " 'splitting" ' operator to [22 



d 2 <$> k 



rk 

-k / dk'<$> k TTZ^Ze i{ - k - h ' )z ^ z + ikN$ k 
Jo 



-kY* I dk'e^e^'^r* + ikN V e Mr " 



(24) 



This term requires careful manipulation: 



Uk = -kj^ f dk'e^'^e^-^r^ + ikNY^e^ 

= -kj^ f" dk'e^e^-^r] -k 2 ^ r^ kr " + ikN ^ e lkr > 

A-J.A JO ' „• 



(25) 



k 2J2r 2 e ikr * +ikNj2e 



ikrj 



3 iv»2 /vi2 

i 3 



-2ik^ 



r 2 e ikr* 
r 2 e ikr] 



ik{N elkr * ~ k2 J2 r * eikr * + ikN e ' 



krj 



ikr? 



2 j. 2 

m 1 3 



hence 



dk 



2tt 



e tkx w k 



-d T 



dy$(y) 



-d x 



x$(x) ( 2 



x-y 
dy$(y) _ d x <£>(x 
x — y $(x) 
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$(x) - (9 x (a;$(a;)) 



(26) 



As it has been shown to be consistent [22], the d x <&(x)/<&(x) term can be 
neglected and we have finally 

«, = -«. [*•(*) f2/*5M)l (27) 
V J x-y J _ 

In the density variable descripton (22], the Jacobian J associated with the 
change of variables to the invariant operators satisfies: 

dX ' n -'^ +W * + I dX '§§)) J = ° 
Since / dx'dn xx >/d®(x') = 0, if follows from and Q2ZD that 



<9$(x) y x — y 

The solution is 

inJ = I dx J dy${x)${y) In \x - y\ ; J = i(r? - r|) 2 = A^(r<) 

in precise agreement with the results of Section 2, e.g., (jHJ). The prefactor in 
(ED is simply the result of the change of variables x = r 2 . 

How does the repulsion amongst the radial eigenvalues express itself as a 
contribution to the potential? This is given by [22] : 



f oo 



loVI uj = / dx x$(x) 



dy $(y) x " 



'o \Jo x y 

Let us introduce a density of radial eigenvalues 0(r) such that 



dx<&(x) f(x) = / 2rdr<b(r 2 ) f(r 2 ) = / dr<p(r) f(r 2 ) 
Jo Jo 

and extend the domain of definition of <f>(r) to the real line by requiring 
d>(—r) = Mr). Then 



'It is associated with a i — j regularization and is, in any case, subleading in l/N. 
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Remarkably, as is the case in the collective field descripton of the singlet 
sector of the single hermitean matrix a local cubic potential is generated 
in the bosonized radially symmetric sector of the 2 + 1 fermions. 



6 More complex matrices 

As is well known, the case of larger number of matrices is of great importance 
in the context of the AdS/CFT corespondence pE], [5], [6]. Of particular 
importance is the case of 3 complex matrices, which are associated with the 
3 Higgs of the bosonic sector of N = 4 SYM. 

Let us consider in general m complex matrices , A — 1, m. Then 

Y, Z A Z A 

A 

is an Hermitean, positive definite matrix. As in the previous subsections, we 
will denote its eigenvalues by rf. The corresponding density variables are: 



$ fc = Tre ifc ^ z b z b = e ikr " ; $(x) = j — e ~ lkx ^ k = S(x - rj) 

i i 

(29) 



One has 



3 In this case the cubic potential is the Thomas-Fermi density term of 1 dimensional 
fermions 
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-f* = ik (z xe «t*).. = a (^»l».zi) (so) 

d(ZA)ij V 0\^A)ij v 'J* 

For f2, the result is identical to the that of the previous section: 

°* = J^t £\ = -"'^^' lt+t ' )4Z - OH 
d{Z A )] j o{^A)ji 

Qxx ' = l^l ^ e ~ ikXe ~ ik ' X ^' = d A' [x*(x)6(x - x')} 
For u>, one obtains 



d 2 ^ 



-k / dk'<$>kRZ ] Ze i{k - k ' )z ^ z + ik mN $ A 
Jo 



d{Z A %d{Z A ) 3l 
= -kJ2 I dk'e ik ' r h i{ - k - k '>h 2 + ik mN ^ e 

ij i 

As described in the previous section, this yields 



(32) 



w x = —d x x$(x) 
and hence the Jacobian satisfies 

d 



dy$(y) | N(m-1) 
x-y 



x 



(33) 



d$(ar) 



In J 



_ 2 Jdy$i 



dy<f>(y)_ N(m - 1) 

y x 



A full treatment of the ensuing Jacobian is beyond the scope of this article, 
but the first term on the right hand of side shows it to be multiplied by 
the Van der Monde type determinant Aj^ R , associated with the radial inter- 
eigenvalue repulsion. 



7 Comments and Conclusion 

By introducing matrix valued polar coordinates, we were able to obtain the 
form of a two hermitean matrix hamiltonian in terms of the eigenvalues of the 
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"radial matrix" (or the eigenvalues of ZZ^) and of a single "angular" unitary 
matrix (Q). For potentials which only depend on the radial eigenvalues, the 
radial Schroedinger equation can be mapped to the Schroedinger equation 
of a system of non-interacting 2 + 1 dimensional "'radial fermions'". We 
obtained the form of the additional potential resulting from the repulsive 
interaction amongst the radial eigenvalues in the density description of the 
system. The presence of repulsive interactions amongst the radial eigenvalues 
has been identified for an arbitrary number of complex matrices. 

The emergence of higher dimensional laplacians from matrices and the 
appearance, as a result of quantum effects, of repulsive inter-eigenvalue in- 
teractions, is potentially of great interest to the AdS/CFT correspondence 
0- jS] and to the general issue of the emergence of gravity (e.g., [30] ) 
and it deserves further study. 

As a note of clarification, clearly the use of the term "fermionization" 
is appropriate to the one dimensional radial wave funtions described in this 
communication. We do not use it in the sense of a full 2 + 1 dimensional 
fermionization. However, even in the presence of " angular" excitations, the 
redefinition ffTSl) would still map the laplacian ([TBI) to the much simpler form 
of (l20j) . added by the same angular contribution. 
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